Vortex line in a neutral finite-temperature superfluid Fermi gas 
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The structure of an isolated vortex in a dilute two-component neutral superfluid Fermi gas is 
studied within the context of self-consistent Bogoliubov-de Gennes theory. Various thermodynamic 
properties are calculated and the shift in the critical temperature due to the presence of the vortex 
is analyzed. The gapless excitations inside the vortex core are studied and a scheme to detect these 
states and thus the presence of the vortex is examined. The numerical results are compared with 
various analytical expressions when appropriate. 
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I. INTRODUCTION 

The achievement of Fermi degeneracy in a confined gas 
of alkali atoms 0, 0, 0, 0j IE IS nas spurred great inter- 
est both theoretically and experimentally in cold atomic 
gases with Fermi statistics. The atomic interactions are 
well-understood and often may be tailored through the 
physics of Feshbach resonances by the application of ex- 
ternal magnetic fields 0, 0, . When the atom-atom 
interaction is attractive, the ground state of a two com- 
ponent gas is predicted to be superfluid at low temper- 
atures 11]. Such a superfluid would provide a unique 
test bed for the study and interpretation of analogous 
but much more complex systems, such as superfluid 3 He, 
unconventional superconductors, and neutron stars. 

One important issue facing the cold atom community 
has been how one would go about actually detecting the 
presence of superfluidity in these systems. Superfluidity 
in Bose-Einstein condensates (BECs) can be inferred ei- 
ther by probing directly the momentum distribution of 
the cloud, the collective modes (where the spectrum is 
strongly shifted relative to the normal phase), or by gen- 
erating quantized vortices (an unambiguous signature of 
the breakdown of irrotational flow) and simply viewing 
the associated "holes" in the particle density 0, Il3| . 
Likewise, for superfluid Fermi gases, the presence of su- 
perfluidity has been shown to give many observable ef- 
fects on the mode spectrum of the gas 0, 0. For 
fermions in the weak-coupling limit, the presence of a vor- 
tex would be very difficult to image directly by looking at 
the density profile, as there is very little depletion of the 
density in the vortex core |16j . However, the quantization 
of angular momentum which is a striking macroscopic 
effect of superfluidity can, as for bosons, be measured 
through the energy shift of the quadrupole modes jl7j . 

Experimental techniques currently limit the tempera- 
ture of trapped Fermi gases to not much less than one 
tenth of the Fermi degeneracy temperature Tp. The su- 
perfluid transition temperature T c of a conventional uni- 
form Bardeen, Cooper, Schrieffer (BCS) superconductor, 



however, is typically lower: T C /T F ~ 0.28e~ Tr/2feF|a| < 1, 
with kp the momentum at the Fermi surface, a the s- 
wave scattering length for low-energy two-body collisions, 
and kp\a\ -C 1 in the weak-coupling approximation where 
BCS theory is valid. A number of schemes to raise T c to 
a value closer to temperatures already accessible with 
dilute Fermi gases have recently been proposed. One 
of these, referred to in the literature as "resonance su- 
perfluidity" involves tuning the scattering length to an 
extremely large value at a Feshbach resonance 
recent experimental results (see for example [53, El l22j |) 
show significant progress using this approach, culminat- 
ing in the production of a Bose-Einstein condensate of 
molcules [23, 13 HH- Another proposal involves load- 
ing the cold Fermi gas into a three-dimensional optical 
lattice [2(j : if the lattice is made sufficiently deep, the 
lowest-lying band will flatten to the point where all of 
the atoms participate in the pairing, as opposed to regu- 
lar BCS theory, where only the small fraction of particles 
close to the Fermi surface are available for pairing. Of 
course, the lattice depth cannot be so great that coher- 
ence across the sample is destroyed, as has been observed 
for bosons in optical lattices [23, Ha H3 ■ The inability 
to experimentally attain very low temperatures in dilute 
gases is probably not fundamental, however. With an eye 
on future experiments, it seems reasonable to explore the 
predictions of a weak-coupling theory of Fermi superflu- 
idity. 

In the present manuscript, we examine in detail several 
properties of the vortex phase of a neutral Fermi liquid 
using a microscopic weak coupling theory. The theoret- 
ical framework is briefly discussed in Section [HJ and we 
present in Section IIIII the details of our numerical pro- 
cedure. Section Hvl is devoted to the calculation various 
thermodynamic quantities of the vortex phase, which are 
compared with the corresponding quantities in both the 
normal state and the superfluid with no vortex. Further- 
more, we demonstrate that the vortex causes a shift of 
the superfluid transition temperature. Finally, in Sec- 
tion we propose a way of observing the vortex through 
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"laser probing" of the quasi-particle states trapped inside 
the vortex core. 



II. THEORETICAL BACKGROUND 

We consider a two component Fermi gas consisting of 
particles with internal quantum numbers a =T>I an d 
mass m a confined in a cylinder of length L and radius 
R. For atomic gases at low temperatures and realistic 
densities, the interactions far from Feshbach resonances 
are characterized by the low energy parameter a which is 
the s-wave scattering length appropriate for the scatter- 
ing between the two specific internal states of the atoms. 
Therefore, only Fermi particles in different internal states 
are able to interact. In our calculations, we assume an 
equal population of the two components iVj = N± so that 
their densities n a are equal. The superfluid phase of the 
gas for a < can be described within mean field theory 
by the Bogoliubov-de Gennes (BdG) equations [3(j 
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• (1) 



Here H KF = - + V{r) + gn a (r) with the low en- 
ergy effective coupling constant given by g = 47rh 2 a/m a - 
The particle density and pairing field are defined as 
n CT (r) = (tpl (r)^ CT (r)) and A(r) = -^(^(r)^^)), re- 
spectively, where V4-( r ) is the usual fermionic field opera- 
tor creating a particle in the internal state a at position r. 
The Bogoliubov wave functions u v (r) and v„ (r) describe 
quasi-particle excitations with energy E v > 0. The ultra- 
violet divergence in the definition of the superfluid gap is 
regularized using the pseudopotential method [3lJ • Since 
the system is essentially homogeneous, the spectrum is 
continuous, so we can use a semi-classical version of this 
scheme also described in [{^ . We augment it to incor- 
porate the effect of the Hartree mean-field gn a in order 
to achieve faster convergence of the solution. The same 
method was used by Grasso and Urban |33| , who present 
a detailed analysis of the convergence properties. 



A. Vortex phase 

The superfluid order parameter is a complex number 
and can thus be written as a real amplitude times a phase 



A(r) = |A(r)|e 



(2) 



The superfluid velocity is then given by the spatial vari- 
ation of the phase of the order parameter 



2m a 



V0(r), 



(3) 



where 2m a is the mass of a Cooper pair. For rotational 
currents, the order parameter must vanish at the center 
of rotation where the superfluid velocity diverges. Far 



away from the core, the flow velocity decreases with the 
distance p from the vortex line as 



2m a p 



(4) 



Here k is the strength of the vortex line. This form of the 
velocity field implies the existence of a region close to the 
vortex axis where the kinetic energy is large enough to 
break the Cooper pairs. Hence the order parameter will 
be suppressed in the vortex core and will heal to its bulk 
value over a length scale governed by the coherence length 
£,bcs{T) = 7ii>F/7rAo(T), with vf = hkp/m a the Fermi 
velocity and Ao(T) the temperature dependent value of 
the bulk gap away from the vortex core [34[ . 

Due to the single-valuedness of the order parameter 
the phase 6 must return to the same value modulo 2tt 
when going around the vortex line. Hence the circulation 
<fv s -d£ is restricted to integer multiples of h/2m a . In 
the present work we will concentrate on vortices of unit 
circulation § v s • d£ = h/2m a . 

In summary, a vortex line represents a topological de- 
fect in the superfluid order parameter, around which the 
superfluid velocity field v s is tangential. The quantiza- 
tion of the circulation represents one the hallmarks of a 
superfluid, and therefore the production and subsequent 
detection of quantized vortices in an ultra-cold atomic 
Fermi gas would be a clear signature for superfluidity in 
the system. 



III. COMPUTATIONAL METHODS 

For a gas confined in a cylinder of radius R and 
length L it is natural to work in cylindrical coordi- 
nates (p,z,tp), where p measures the perpendicular dis- 
tance from the symmetry axis, z is the axial coordinate, 
and if is the azimuthal angle around z. In this coor- 
dinate system the order parameter can be written as 
A(r) = |A(p, z)\ exp(— imp), with k = corresponding 
a phase with no vortex, and k = 1 for a singly quan- 
tized vortex along the axis of symmetry. The mean-field 
density is rotationally invariant: n a (r) = n a (p,z). 

Assuming free motion along the cylinder axis, and im- 
posing periodic boundary conditions at z = ±L/2, we 
write for the quasi-particle modes 

= „-V2 



■^(r) = p I U nm kAp) 



2tt VZ 



fe, (P) 



(5) 



The allowed values of the angular momentum quantum 
number are {m = 0, ±1, ±2, . . .}, and k z — 2tt£/L, with 
{£ = 0,±I,±2, . . .}. The radial functions (u nmkz ,v nmkz ) 
are taken to be real. With these definitions the BdG 
equations Q become 



A(p) 

-H. 



ikAp) 



= E n 



U„mkz (P) 
Vnmk z (P) 



,(6) 
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where 



Hm 



2m a 



cP_ 
dp 2 



1/4) 



+ gn u {p) -/z. 



(7) 

These are the equations we solve self-consistently through 
an iterative procedure. 

By exploiting the symmetry of the BdG equations Q , 
we can identify a negative energy solution with angular 
momentum m with a 'positive energy solution with an- 
gular momentum — to — k. We can therefore generate 
the entire positive energy spectrum by solving Eq. © for 
?7i > only, and using the transformation 



E„ 



-E, 



v 

-ul 



(8) 



to find the eigenstates with to < 0. 



A. Discrete Variable Representation 

The BdG equations in general must be solved numer- 
ically. Some of the effects of the vortex that we are in- 
terested in, such as the associated shifts in the critical 
temperature T c and in the ground state energy of the 
gas, are quite hard to calculate numerically as they are 
very small compared with the corresponding bulk val- 
ues. For example, to obtain the vortex energy one needs 
to subtract two large numbers (the ground state energy 
of the gas with and without a vortex) to get a small 
number. This requires a very accurate numerical scheme 
to solve the BdG eqns. Such a scheme is provided by 
the Discrete Variable Representation (DVR) which re- 
cently enabled the microscopic calculation of the vortex 
energy ^6] . DVRs are representations on a basis of func- 
tions localized about discrete values of the coordinate. 
This renders local functions of the coordinate operator 
approximately diagonal within the DVR basis, making 
DVRs ideally suited for solving self-consistent problems 
like the present one, where the matrix elements of the 
pairing and Hartree fields (local functions) have to be 
evaluated at each iteration. In addition the representa- 
tion of the kinetic energy operator is exact. The litera- 
ture on DVRs is extensive and we shall only convey the 
central points here. A detailed review of the framework 
can be found in [3^, |3|J . 

A DVR exists when there is both a spectral basis of 
N functions, (f> n (x), orthonormal over an interval [a,b] 
with weight function w(x) and a quadrature rule with iV 
points x k and weights w k 

,6 N 
(f\g)= dxw{x)f(x)g{x) = ^w k f(x k )g{x k ). (9) 
Ja fc=i 

This enables a set of coordinate eigenfunctions 
{tpi(x), i = 1, N} to be defined with the property 



Iw(Xi) 



(10) 



We expand the unknown functions ipi (x) on the basis 



N 



Ipiix) = ^2 <t>n(x){<l>n\lpi), 



(11) 



and use the quadrature rule I© and (jl()|l to evaluate the 
expansion coefficients. The coordinate eigenfunctions are 
then given by 



JV 



jJijx) = y / w(x)w i (t) n (x)(f> n (x i ). (12) 

n=l 

Since the ipi(x) diagonalizes the coordinate operator, the 
matrix element of any operator 0(x), which is a local 
function of x, is approximately diagonal within the DVR 



(13) 



the approximation being due to the use of a truncated ba- 
sis. Furthermore, since the DVR involves an underlying 
spectral representation, it is possible to evaluate matrix 
elements of parts of the Hamiltonian exactly, if the <j> n {x) 
are chosen to be the eigenfunctions of the corresponding 
operator (for example the kinetic energy). 

For the problem of quantization in a cylinder the cylin- 
drical Bessel functions form an ideally suited basis for the 
DVR as suggested in reference [37j • They are orthogonal 
over the range [0, R] 



R 5. 
dp p J m (ki m p)J m {kj m p) 



(14) 



where the coordinate normalization constant is given 

by m 



(15) 



Similarly, the Bessel functions are also orthogonal in mo- 
mentum space: 



K g 

dk k J m {kpim)Jm{kPjm) — • 
Wi m 



with the momentum normalization 

2 



(16) 



(17) 



The spatial and momentum grids are pi m — Zi m /K, and 
hm = Zi m /R, respectively, where {z irni i = l,...,N} 
are the zeros of the Besselfunction of order to, defined 
through J m {zim) = 0. This is a consequence of the 
boundary condition which states that the wave function 
must vanish at p — R. Note that since fc/vm = z^ m /R = 
K mi and pn = ZN m /K m = R, the maximum momen- 
tum and the maximum value of p are not independent, 
but are inversely related to each other by the relation 
RK m = ZN m . It was shown in that a quadrature 



4 



rule can be associated with these grid points, provided 
weights are chosen to be Wi m (w' im ) for integration over 
the spatial (momentum) variable. In general there will 
be one spatial and one momentum grid associated with 
each value of the angular momentum m. 

With the Bessel function quadrature in place we can go 
ahead and construct a DVR basis. As our orthonormal 
basis functions we choose 



350.8 



MP) = VKmJrnihmP), 



(18) 



where the yj w' im is necessary to ensure that the basis set 
is orthonormal, i.e. (4>k\4>i) = $kl- From (|12|l we thus 
have for the coordinate eigenfunctions 

im^nm (^Virn p) Jrn, \^nm Pirn) • 

(19) 



The radial functions (u r . 



he 



panded in terms of the coordinate eigenfunctions, i.e. 
u nm k z (p) = J2i limi>im{p)- The BdG equations will then 
be a set of non-linear equations for the expansion coeffi- 
cients 7i. Due to the properties of the coordinate eigen- 
function the value of the radial function on the grid points 
is simply u nmk __(Pi m ) = Ji m y/pim/w im . 

We conclude this section with two important remarks. 
While the transformation from the spectral basis to the 
coordinate eigenfunctions is not strictly unitary, the nu- 
merical procedure is nonetheless well defined, as the 
transformation can be made unitary in the limit of large 
N [37|]. Secondly, although it appears that a separate 
grid is needed for each m value, we have found that in 
practice only two grids are needed, one based on Jo for 
even m and one based on J\ for odd m. Since u v and 
v n for a vortex state correspond to wavefunctions which 
differ by one unit of angular momentum, they will be 
represented on different spatial grids. Fortunately, inter- 
polation is trivial in the DVR method. To interpolate 
from the m = to the m = 1 grid amounts to multi- 
plying the vector of expansion coefficients 70 with the 
transformation matrix given by Bij = ipio(Pji)- The re- 
verse transfomation is b\- — 4>n{Pjo)- For the purpose 
of solving the BdG equations the mean-fields are only 
represented on the odd m grid. 



IV. THERMODYNAMICS 

In this section, we present results for various thermo- 
dynamic quantities of the vortex phase obtained by solv- 
ing the BdG-eqns. numerically as described above. All 
calculations were done for a fixed N a = 28000. The ra- 
dius an d length of the box were taken to be 28.5 ^m 
and 11.4 /j,m, respectively. For 6 Li the scattering length 
is —2160 ao, which gives a bulk value of the transition 
temperature T c q = 0.01 fiK, and a Fermi temperature of 
Tp = 0.70 fiK for the chosen density. 
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FIG. 1: Plot of the free energy per unit volume of N a = 
28000 fermions in the normal, and superfluid phase with and 
without a vortex as a function of temperature. The solid 
and dashed vertical lines represent analytic expressions for 
the vortex and condensation energies, respectively. 

In figure \I\we plot the free energy (H) — TS as a func- 
tion of the temperature T. The entropy is found as 

S = -k B ^ [f(E v )\nf(E v ) 

+ (1 -/(£„)) ln(l -/(£,,))], (20) 

since the quasi-particles in our mean-field approach form 
an ensemble of non-interacting fermions 30]. We have 
calculated the free energy for the vortex phase, the 
superfluid phase without a vortex, and for the nor- 
mal phase. All have been normalized to 7T C 2 , where 
7 = 2n 2 N(Q)k 2 B /3, and N(0) = 3n a /2e F is the den- 
sity of states per unit volume (for a single component) 
at the Fermi energy in the normal phase [3!j. For 
T = 0, the condensation energy density of the super- 
fluid without a vortex with respect to the normal phase 
is E cond /V = -JV(0)Ag/2, with A = % e -z e -*/^ F \a\ 
the bulk value of the superfluid gap. This condensation 
energy is indicated on the figure and we see that there 
is good agreement with the numerical results. Further- 
more, the vortex energy per unit axial length for T = 
due to the loss of condensation energy in the vortex core 
and the kinetic energy of the supercurrent around the 
core is |l7j 



irh 2 n a 
2m a 



In 



D- 



R 



£bcs(0) 



(21) 



The constant D was determined numerically in ref. [l6| 
to be D w 2.5. This expression for the vortex energy is 
also indicated on the figure. 

As can be seen from Fig. |3 the critical temperature 
for the vortex phase T cv is lower than that of the bulk 
superfluid phase without a vortex T c q. For the specific 
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FIG. 2: Chemical potential in the normal, and superfluid 
phase with and without a vortex, as determined by the con- 
straint that N a = 28000. We attribute the low temperature 
behavior of the normal phase chemical potential to shell ef- 
fects due to the finite volume |40j. For the vortex state the 
transition temperature is shifted downwards. 



parameters used, the difference is 1 — T cv /T c o ~ 0.1. 
This difference can be understood as follows: The vor- 
tex phase becomes unstable with respect to the normal 
phase when the extent of the vortex core becomes com- 
parable to the radius R of the system. Since the size of 
the vortex is 0(£bcs)> we can estimate T cv from the con- 
dition £ B cai2>) ~ OCR). Using A (T) a 1.7A (0)(1 
t/t U/2 Tori 



T/Tcq) 1 / 2 |30j for < 1 - T/T cQ < 1, this yields 



6Tc 



1 



£bcs(0) 5 



(0 



R 2 



-Oil 



(22) 



where ai is a number of order one. We now test this 
expression and determine the constant ct\ by numerically 
calculating the shift in the critical temperature ST c /T c q 
due to the presence of a vortex for various radii of the 
system. The result is shown in Fig. [31 We find that we get 
reasonable agreement with Eq. (|22[) as £bcs /i? — > with 
a coefficient oti ~ 2.3. So one can understand decrease in 
T c due to the presence of the vortex as a finite size effect 
which scales as £bcs(0) 2 /-R 2 - 

In Fig. 0| we plot the temperature dependence of the 
heat capacity cy = (T/V)dS/dT per unit volume of the 
system. Again, we show for comparison results both for 
the system in the normal phase, in the superfluid phase 
without a vortex, and in the vortex phase. For a two com- 
ponent gas in the normal phase, we have cy normal = "fT 
for T — > 0. The heat capacity for the superfluid phase 
without the vortex is exponentially damped by a factor 
exp(— j3Ao) for T«T C due to the gap in the energy spec- 
trum [33 ■ Figure0]on the other hand shows that the heat 
capacity in the vortex phase cy vortex depends linearly on 
T for low temperatures. This linear T-dependence is due 
to the presence of so-called core bound states in the vor- 
tex phase. These are single-particle excitations which 
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FIG. 3: The shift in the superfluid transition temperature for 
the vortex state relative to a bulk superfluid with no vortex, 
as a function of the radius of the confining cylinder at fixed 
density. 



are spatially localized in the vortex core where the gap is 
small. The energy of the core states is in general less than 
the bulk gap energy Ao and they exist only for angular 
momentum quantum numbers m > [30l |. This corre- 
sponds to a quasi-particle current around the vortex core 
in the opposite direction to that of the vortex current. In 
a detailed analysis it was found that the energy spectrum 
of the lowest bound core states with < m <C ^f^bcs 
for T — is essentially gapless and given by 



E mkz ~ (m + 1/2) 



ep sin 8 



(23) 



where k~ = Uf cos 9 and h{9) is a function of order 
unity In Fig. [S] we plot the lowest quasi-particle 

energies as a function of m for k z = obtained from a 
numerical solution of Eq. JJ]J. The gapless branch asso- 
ciated with the core states with energies less than Ao is 
clearly visible. The T — density of vortex states per 
unit volume is calculated by integrating Eq. I|23|l over k z 
which yields 



N v (e) 



(24) 



for < e < A where a 2 ~ 0(1) g^]- Thus, the den- 
sity of bound core states per unit volume is the same, 
apart from a factor CV2, as that of a cylindrical region of 
a single component gas in the normal phase with radius 
£bcs and length L. From this we conclude that the low T 
heat capacity per unit volume of the gas in vortex phase 
cy vortex associated with the core states is 



Cy vortex ^ CVnormal^2 



£bcs 

H 2 



(25) 



explaining the linear T-dependence of cy vortex observed 
in Fig. Q] A fit to the numerical data yields a.2 ~ 2. 



6 



2.5 



1.5 



0.5 



0.03 




0.05 0.1 0.15 0.2 



q! » « » » * 



0.2 0.4 0.6 0.: 

T/T „ 



1.2 



FIG. 4: Plot of the specific heat per unit volume in the nor- 
mal, and superfluid phase with and without a vortex. The 
inset shows the low temperature behavior for the vortex state 
and the superfluid state without a vortex (same symbols as 
in Figs. □ and 0. 
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FIG. 5: Energy spectrum for the lowest quasi-particle states 
in a superfluid with a vortex (•) and the vortex free state (x) 
at T — 0. For clarity only the energies of states with k z = 
have been plotted. There are branches of bound states for 
several values of k z . 



We remark that a linear contribution to the heat capac- 
ity has been observed for a superconductor in the mixed 
state IH. 



V. LASER PROBING OF THE VORTEX PHASE 

Vortices are now routinely created in dilute BECs 
where they can be detected by direct imaging of the 
cores, in which the density is significantly suppressed. 



Unfortunately, such a procedure would be very difficult 
to implement successfully for a dilute superfluid Fermi 
gas, where there is no significant depletion of density in 
the vortex core One way to observe the vortex is 

to measure the shift in the quadrupole mode frequencies 
which is directly proportional to the angular momentum 
per particle h/2 associated with the supercurrent around 
the vortex core |17| . 

In the present section, we investigate the feasibil- 
ity of detecting the bound quasi-particle states in the 
vortex core through a recently proposed laser probing 
scheme |44l |45| . The laser probing scheme is similar to 
scanning tunneling microscopy (STM) on a supercon- 
ductor in that it relies on induced tunneling between 
a superfluid and a normal phase 0, ^?|]- Whereas a 
STM probe uses a bias voltage to transfer population 
across a superconducting-normal interface existing be- 
tween the normal microscope tip and the superconduct- 
ing substrate, the laser probe instead creates an effec- 
tive interface by coupling different internal states of the 
atoms by laser fields. Specifically, a spin state | |), which 
is Cooper paired with the state | J.) is coupled via laser 
field to a third state \e) that has been chosen such that it 
does not participate in the pairing (either it does not have 
strong attractive interactions with the two other states or 
the disparity in chemical potentials is too large). Hence, 
the |e) atoms define the normal part of the interface. If 
the detuning of the laser from the atomic transition is 
5 = iu a — wl, where ujl is the laser frequency and to a 
the frequency splitting between the level | f) and |e), the 
rate of change in the population of the |e) state (tunnel- 
ing current) I = —(N e ) is |45| 



1(5) 



2tt » - 



^!J(r)«,(r)$;(r) 



x {f(E v ) - f(E n )]8{E v - e„ - S) 

2 

d 3 r n(r)v;(r)K(r) [I - f{E v ) - /(&)] 

x 6(E v +S n + 8). (26) 

Here S = /i e — /i + 5 = A/x + 5 is the effective de- 
tuning, fi e the chemical potential of the \e) atoms, and 
$ n their single particle wave functions with energy £„; 
fix) = [exp(/3x) + is the Fermi function and f2(r) 
the Rabi frequency. In the present analysis, we assume 
for simplicity that the |e) atoms are non-interacting such 
that their wavefunctions <!>„ are the eigenstates of the 
confining cylindrical box. We consider the case of a con- 
stant laser profile O(r) = f2. This gives the selection rule 
kj = k e where kf is the momentum of an | f ) atom cou- 
pled by the laser beam to an |e) atom with momentum 

Let us now consider how the laser probing method can 
be used to probe the presence of the core states. We 
examine two opposite cases of interest: The case when 
there are initially no |e) atoms present (N e = 0) and the 
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case where there initially are an equal number of | f ) and 
|e) atoms present (JVf = N e ). 

From Eq. ({26(1 it is straightforward to show that for the 
total current we have J dSI(5) oc N e — Np That is, the 
net current from the |e) atoms to the | f) atoms is pro- 
portional to the difference of initial populations between 
the two hyperfine states. Likewise, the total current from 
the core states trapped inside the vortex is clearly pro- 
portional to the total number of core states N cs . Thus, 
when there initially are no \e) atoms present (N e = 0) 
the spectral weight of the current due the core states as 
compared to the total current observed scales as N cs /Np 
Using N cs ~ N v A a nR 2 L with N v given by Eq. lf2"I)l. one 
obtains that the current from the core states divided by 
the total current scales as Aoe^ 1 ^B CS -R -2 <C 1. Thus, the 
signal from the core states is completely overwhelmed by 
the bulk signal coming from the current out of the whole 
Fermi sea of | |) atoms. We therefore conclude that it is 
most likely not possible to probe the core states starting 
with initially no |e) present. This conclusion is supported 
by numerical simulations. 

Let us therefore consider the case when there ini- 
tially are an equal number of | f) and |e) atoms present 
(iVf = N e ). In that way, the bulk signal of transitions of 
| T) atoms deep within the Fermi sea is Pauli blocked due 
to the presence of the |e) atoms since we have the selec- 
tion rule kj = k e . One can then show from Eq. 1|2()|) that 
the total signal scales as J d5\I(5)\ oc N^Ao/ep, i.e. the 
current is proportional to the total number of Cooper 
pairs. Thus, the bulk signal is suppressed by a factor 
Aq/cf compared to the case when there are no |e) atoms 
present simply due to the Fermi blocking effect. The 
current due to the vortex core states should therefore be 
easier to observe as it is not overwhelmed by a huge back- 
ground signal. In Fig. [5] we plot the T = laser probing 
current 1(6) for the case when N e = Np The effect of 
the Hartree field gn a is primarily to shift the entire pro- 
file to lower detunings 5 since it shifts the energies of 
the | |) atoms by the amount gn a whereas the |e) atoms 
are assumed non-interacting. In the plot we have explic- 
itly eliminated this overall shift for reasons of clarity. We 
plot the current both when there is no vortex present and 
when there is a vortex. In the case of no vortex present, 
the current given by Eq. I|26|l at zero temperature can be 
shown to be 

I = ±nn 2 p(5)e(S*-At)^, (27) 

where ± corresponds to 5 > and 5 < respectively 
and p(S) = V{Al/5 -5 + 2p) 1 / 2 /2n 2 0,111. From Eq. 
(|2*7|l it follows that there is no current for detunings with 
— A < S < A . This can be interpreted as the laser 
signal has to provide a minimum energy Ao to break 
a Cooper pair and generate a current. Equation (|27|l is 
also shown on Fig.[f)]and we see good agreement with the 
numerical result when there is no vortex present. Note 
that since the numerical calculations use a Lorentzian of 



width r = O.OlAo instead of S(x) functions in Eq. 12611 . 
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FIG. 6: The tunneling current as a function of detuning (in 
units of the bulk value of the gap) for tunneling into filled 
state from both a vortex state and a superfluid without a 
vortex. For comparison Eq. I|27|l is also plotted. The profiles 
have been shifted to compensate for the Hartree mean field 
shift gricr of the energies of the | T) atoms. If the | j) atoms 
are in the normal state no current flows due to Pauli blocking. 



we have convoluted Eq. (|27|l accordingly. We see that 
the signal when there is a vortex present is markedly 
different from the case with no vortex. In particular, 
there is a significant current for \5\ < Aq. This current 
is directly due to the presence of the core states which 
have a pairing energy less than Ao- The signal from the 
vortex phase is finite for 5 ~ reflecting the fact that the 
energy spectrum of the core states approximately given 
by Eq. (|2~5|l is essentially gapless. Thus, the existence of 
core states bound in the vortex is reflected in the current 
profile 1(6). 



VI. CONCLUSIONS 

We have studied the properties of a single vortex in 
a neutral superfluid with Fermi statistics using a micro- 
scopic weak coupling theory. The effect of the vortex on 
the free energy and the heat capacity of the system was 
examined and we provided various analytical expressions 
which agrees well with the numerical results. The vor- 
tex gives rise to the presence of core states bound in the 
vortex core. We examined the spectrum of these states 
and also suggested a way to experimentally detect them. 
Apart from being of interest theoretically, it is not un- 
likely that our results will have experimental relevance in 
the near future due to the recent impressive experimental 
progress within the field of atomic Fermi gases. 
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